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We use the tunneling formalism to calculate the Hawking radiation of massive particles. 
For E > m, we recover the traditional result, identical to the massless case. But E < m par- 
ticles can also tunnel across the horizon in a Hawking process [23| . We study the probability 



for detecting such E < m particles as a function of the distance from the horizon and the 
energy of the particle in the tunneling formalism. We derive a general formula and obtain 
simple approximations in the near-horizon limit and in the limit of large radii. 

I. INTRODUCTION. 

n 

Hawking radiation [1] consists in the emission of pairs of quanta from a black-hole horizon: one 
with "positive energy" (positive co-moving frequency) towards the exterior and one with "negative 
energy" (negative co- moving frequency) towards the interior of the black hole. The emission has a 
thermal spectrum as seen by an asymptotic observer at infinity, with a temperature Th determined 
by the derivative of the free-fall velocity u at the horizon r/^: Th = h\^\^_^^ /2tt. This result can be 
obtained very simply in the semiclassical method. This tunneling description of Hawking radiation, 
first introduced by Volovik (see also 31 for the Parikh-Wilczek version, {4] for its extension to 
E > m massive particles, and the recent 5|] for charged E > m massive particles), indicates that 
the Hawking process can be understood as the quantum tunneling across the horizon between 
classical trajectories on both sides of the horizon. The semiclassical method was also applied to 
Hawking radiation from rotating [6] and charged j?! black holes, and to several higher-dimensional 



and other more exotic black hole geometries (see e.g. [s-llll] for some recent examples), as well as 



to related phenomena such as the Zel'dovich-Starobinsky effect 13] (see [3]), and the Unruh 
effect lisl (see e.g. Q). 

We will focus on massive particles in a Schwarzschild geometry. If one is interested only in the 
spectrum as measured at infinity, then it is straightforward to arrive at the usual conclusions jl?! . 
namely that there is a threshold E > m for massive quanta, but otherwise, the spectrum does not 
differ from that of the standard massless case. 
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However, although no E < m particles arrive at infinity, they are nevertheless radiated from 
the horizon, see also [Q]. We will calculate the probability for detecting such E < m particles as a 
function of the distance from the horizon and the energy of the particle. We largely follow Volovik's 
derivation and notation Q] (see also 13 ) , focus on purely radial movement in a Schwarzschild 
geometry and set ?i = c = 1. 

II. GENERAL CALCULATION. 

We start from the Painleve-Gullstrand-Lemaitre (PGL) form of the Schwarzschild line element: 
ds"^ = g^ydx^'dx'' = -(1 - v'^{r))dt^ - 2v{r)drdt + dr'^ (1) 



where v{r) = —^/r^/r is the free-fall velocity and the horizon or Schwarzschild radius. The PGL 
metric is stationary and moreover, unlike the Schwarzschild form, it is regular across the horizon. 
This makes it particularly suited for the tunneling description of Hawking radiation. 
For a massive particle, one has 

9^y^^Pu = -vn^ (2) 
with = {—E,pi), leading to the energy-momentum dispersion relation 

m^+p^ = (^_p. v)2. (3) 
Since v represents the free-fall velocity, it is natural to write this as 

S = So + P • V, (4) 
and interpret E as the Doppler-shifted energy (the energy in the "black hole rest frame"), which 



is a conserved quantity, and Eq = \J m? + as the energy in a co-moving reference frame. This 
interpretation is reinforced by the observation that the PGL metric is also the natural metric in 
analogue gravity pi], where it describes the propagation of sound waves or other perturbations on 
a background fluid moving with a velocity v. 

We want to calculate the influence of the mass on the tunneling rate of particles from the black 
hole up to a detector (where the particle can be captured in a bound state, see also at an 
arbitrary fixed radial distance R. The tunneling probability W is determined by the imaginary 
part of the action S along the semiclassical trajectory: 

W oc exp[-2Im5], (5) 
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where 

ImS = Im J pr{r)dr (6) 

and Pr{r) follows from the dispersion relation The final result will be of the form 

W{E) oc exp[-2Im5i] exp[-2ImS2] (7) 

with ImSi a contribution for tunneling through the horizon, and ImS'2 an additional action in case 
there is a second classically prohibited region beyond the horizon. We certainly expect ImS'2 7^ 
when E < m, since the massive particle is then classically forbidden [see Q] in flat spacetime 
[v{r — )• 00) = 0]. 

From ([3]), we obtain 

P = + -r^Vm^v^ -1) + E^ (8) 

1 — 1 — 

= Pl+P2 (9) 

Note that, in order for the semiclassical formalism to be valid, the action must be large (Im5 S> 
1). In particular, this implies that we consider the case mr^ ^ 1. In the opposite limit mrh <C 1, 



it was found in 



18( 1 that the created bound particles are localized mainly on resonant levels outside 



the black hole with large occupation numbers. 

A. Si: Tunneling through the horizon. 

For pi, we shift the contour of integration to the complex plane and apply the standard residue 
theorem f{z) = 27rzRes/. There is one pole along the radial path: v = —1 at r = r/^, so 

-Ev E , . 

Respi = Res- —, r = — - — r, (10) 

{l + v){l-v) 2v'{rhy ^ ' 

where the prime denotes djdr. We recover the standard Hawking result, as for a massless particle 
to reach infinity: 2Im5'i = E/Th, with Th = |f'(r/i)|/27r. The presence of a mass has absolutely 
no influence on the probability of tunneling across the horizon. 

B. 52 : Tunneling towards the detector at a distance R. 



P2 will only have an imaginary contribution when 

^2 < m^{l - v"^). 



(11) 
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For a Schwarzschild black hole {v = r/j/r), this leads to 



E <E,{R) = m[l-'^\'^ (12) 



R 

for given R. Alternatively, it can be expressed as a condition on R for given E: 

TIT?' 1 

^ > "'^^^^ = ^^123:^'''^ = rr^"" ^^^^ 

For E > Ec (or i? < Tc), there is no additional barrier. Once the particle has tunnelled through 
the horizon, it can freely propagate up to the detector at R. At sufficient distance from the 
black hole {v — )• 0), the condition for a second imaginary contribution reduces to < m, as we 
anticipated. This means that, for E > m, the mass term does not cause any additional tunneling 
factor, independently of R, and we recover the standard result for massless particles all the way 
to r — )• oo. For E < m, however, the particle is created with an energy E which is insufficient to 
escape all the way to oo, so it will encounter a second barrier as v decreases (i.e., as the particle 
moves towards flat spacetime). 

In case of (jlip . we take the (positive) imaginary part and obtain 



R 



1 



dr- —\lvrfi(\ — Tji/r) — E^ (15) 

where the last expression is specific for a Schwarzschild profile. 

Note that the logarithmic divergence for r — )• r/^ (u — )• 1) is avoided because rc > r^- 

Imp2 is shown, together with Im5'2, for various values ol E/m in Fig(TJ For E /m < l/\/2, Imp2 

increases to a maximum /E at r = rhm? / {m? — 2E'^) before decreasing to the asymptotic value 



\/m? — E"^ for r — t- oo. For E/m > 1/^2, the maximum disappears and the increase is monotonic. 



III. LIMIT CASES 

Simple analytic results can be found in the following interesting limit cases. 
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FIG. 1: FiglU Semiclassical barrier Imp2 (thick lines) and tunneling action ImS'2 (dashed lines) beyond the 
horizon for various values of E/m < 1. 




FIG. 2: Figl^l Tunneling action IniiS'2 (thick lines) and approximation (|16p for R ^ rh (dashed lines) as a 
function of E/m. 

A. Limit R:$> rh 

When the detector is very far from the horizon, the integral is dominated by the contributions 
where v ^ and we can immediately write 



pin. 

ImS'2 = \/ — E"^ / dr R\/ m? — E"^ 

Jtc 

where we have assumed rc, see Figj2j The global tunneling probability 



(16) 





'-E' 




oc exp 


[th\ 


exp 



(17) 



decreases exponentially with the distance R to the detector, and the energy difference ^/m? — E'^. 
Note that this result does not depend on the Schwarzschild profile, but is generally valid as long 
as R lies sufficiently far away from the horizon for the integral ()14p to be dominated by the region 
w <C1. 
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FIG. 3: Figini Tunneling action ImS'2 (thick lines) and approximation (|T^ for E <^ Ec (dashed lines) as a 
function of E/m for different positions R of the detector near the horizon. 



B. Limit R 



rh 



Another interesting limit is f — )• 1, i.e. — — 1^1. In this near-horizon limit, there is a non-zero 
contribution only for E <^ m, see (jl3p . There will then be a strong barrier almost as soon as the 
particle crosses the horizon (see the case E/m = 0.1 in Fig{T] above), since it has practically no 
energy to sustain its own mass. 



Starting from (jlSp . we write r /r^ — 1 = <C 1, and obtain 



ImS'2 = Erh 



dy 



(18) 



where we have used <C 1 and 



r ^ rh. 



Two limit cases yield a simple interesting result. 



R 



1 < 1 (i.e. 



R 



1 > 



For given R, this implies E"^ /E1{R) <^ 1. The integral is then dominated by the region where 
y » 1 so 

1/2 



ImSo 



2mrh 1 

Jh 



(19) 



The second tunneling action Im5'2 becomes independent of the energy E in this limit, as illustrated 
in FigH 

l") - 1 « 1, (i.e. ^ - 1 



n ^ R 



For given R, this condition is equivalent to 1 — ^ <^ 1. Now y ~ 1 in the whole integration 
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FIG. 4: Fig|4j Tunneling action ImS'2 (thick lines) and approximation (|22|) for E ^ Ec (dashed lines) as a 
function of E/m for different positions R of the detector near the horizon. 



region, so 



ImS'2 w -r-/, — 7 



2 m" 



^-1 



3/2 



2 ^ 



2 \ 3/2 



R 

-mrh ■ 
3 V^/i 



^2x3/2 



(20) 
(21) 
(22) 



where we have used R ^ r^- 

The prefactor is similar to the previous result, but now there is an additional suppression with 
decreasing E, see Fig SI 



IV. CONCLUSION. 



Massive particles tunnel across the horizon at exactly the same rate as massless particles of the 
same energy E, even HE < m. Such massive particles with E < m do not reach infinity, but 
nevertheless have a non-zero probability of being detected at any finite distance {R — r^) from 
the horizon. For a detector very close to the horizon, or very far from the horizon, the detection 
probability reduces to simple expressions in terms oi E/m and R/r^. Different laboratory systems 
in which this behaviour of massive particles near a black hole could be simulated in the context 
of analogue gravity were suggested in [2^. An interesting prospect for future work is that the 
double barrier could cause a resonance mechanism, thereby leading to peaks of strongly increased 
tunneling probability. Such resonant Hawking radiation in the presence of double-barrier structures 
was recently studied for a BEC-based analogue system in j24l] . 
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